AERONAUTIC  SYMBOLS 

1=  FUNDAMENTAL  AND  DERIVED  UNITS 


Symbol 

Metric 

English 

Unit 

Abbrevia¬ 

tion 

Unit 

Abltroria- 

tion 

Length _ 

meter _  _ 

ft  (nr  mi) 
sec  (or  hr) 
lb 

Time _ 

second. _  _ _ 

Force _ 

F 

weight  of  1  kilogram _ 

kg 

weicht  of  1  pound _ 

Power _ 

P 

horsepower  (metric) _ 

lip 

mph 

Speed _ 

V 

fkilonieter.s  per  hour _ 

knh 

miles  per  hour.  - .... 

tmeters  per  .second _ 

iiip.s 

feet  per  second _ _ 

fps 

w 

Q 

m 

1 


5 

6 
G 
b 

c 

A 

V 

q 


L 


D 

Do 

Di 

O, 

c 


2.  GENERAL  SYMBOLS 

Vreig}it=7n/7 

Standard  acceleration  of  gTavity=9.80GG5  m/s" 
or  32.1740  ft/sec® 


Mass=~ 

g 

Moment  of  inertia=?n/i;^.  (Indicate  axis  of 
radius  of  gyi-ation  Jc  bj’  proper  subscript.) 
Coefficient  of  viscosity 


r  Kinematic  viscosity 

p  Density  (mass  per  unit  volume) 

Standard  density  of  dry  air,  0.12497  kg-m"'-5"  at  IG”  C 
and  7G0  mm;  or  0.002378  Ib-ft"*  see- 
Specific  weight  of  "standard”  air,  1.2355  kgim?  or 
0.07G51  Ib/cu  ft 


3.  AERODYN.^MIC  SYMBOLS 


Area 

Area  of  wing 
Gap 
Span 
Chord 

Aspect  ratio,  -g 
True  air  speed 
Dynamic  pressure,  -^pT^ 


Lift,  absolute  coefficient  (7t=^, 

Drag,  ai)soluto  coefficient 
Profile  drac,  absolute  coefficient 

•'go 

P‘ 

qS 

Pi 

qS 


Induced  drag,  absolute  coefficient  Cd,= 
Parasiie  drag,  absolute  coefficient  Cup 


Cross-wind  force,  absolute  coefficient  Cc~ 


C 


ip,  Anglo  of  setting  of  wings  (relative  to  thrust  line) 
it  Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Q  Resultant  moment 

£1  Resultant  angular  velocity 

E  Reynolds  number,  p^  where  Z  is  a  linear  dimen- 

M  _ 

sion  (c.g.,  for  an  airfoil  of  1.0  ft  chord,  100  mph, 
standard  pressure  at  15°  C,  the  corresponding 
Reynolds  number  is  935,400;  or  for  an  airfoil 
of  1.0  m  chord,  100  mps,  the  corresponding 
Reynolds  number  is  6,865,000) 
a  Angle  of  attack 

f  Angle  of  downwash 

ofo  Angle  of  attack,  infinite  aspect  ratio 

a.  Angle  of  attack,  induced 

Oa  Angle  of  attack,  absolute  (measured  from  zero- 

lift  position) 

•y  Flight-path  angle 
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APPLICATION  OF  THE  METHOD  OF  CHARACTERISTICS  TO  SUPERSONIC  ROTATIONAL  FLOW 

By  Antonio  Kerri 


SUMMARY 

.-1  system  for  calculating  the.  physical  properties  of  svpersonic 
rotational  flow  with  aria!  symmetry  and  supersonic  rotational 
flow  in  a  twoielimensional  field  was  determined  by  use  of  the 
characteristics  method.  The  sy.stem  was  applied  to  the  study 
of  external  and  internal  fiow  for  supersonic  inlets  with  axial 
symmetry.  For  a  circular  conical  inlet  the  shock  that  occurred 
at  the  lip  of  the  inlet  became  .stronger  as  it  approached  the  axis 
of  the  inlet  and  became  a  normal  shock  at  the  axis.  The  region 
in  which  strong  shock  occurred  increased  with  the  increase  of 
the  angle  of  internal  cone  at  the  lip  of  the  inlet.  For  an  inlet 
with  a  central  body  the  method  of  characteristics  was  applied  to 
the  design  of  an  internal-channel  shape  that,  theoretically,  results 
in  very  efficient  recompression  in  the  inlet;  it  was  shown  that 
if  an  ejfuser  is  connected  with  the  diffiuser  a  body  of  revolution 
with  very  small  shock-ware  drag  can  be  determined. 

INTRODUCTION 

TIic  characterist  ics  mctliod  for  the  (Ictcrminatioii  of  super¬ 
sonic  phenomena  was  first  used  hy  Prandtl  and  Busomann 
for  two-dimensional  flow  (references  1  and  2).  For  flow 
with  axial  symmetry  Frankl  (reference  3)  used  the  method 
of  characteristics  for  determining  the  shape  of  a  supersonic 
circular  efliuser  witli  uniform  exit  velocity,  and  Ferrari 
(references  4  and  5)  independently  used  the  charai^teristics 
method  for  determining  supersonic  phenomena  for  every 
type  of  boundary  condition.  Subsequently  Gudei'ley  (refer¬ 
ence.  6)  and  Sauer  (reference  7)  transformed  the  system  pro¬ 
posed  by  Frankl  and  Ferrari  and  obtained  a  different 
analytical  solution  of  the  problem.  In  all  applications  the 
hypothesis  of  potential  flow  was  made;  therefore  the 
equation  of  potential  flow  was  used. 

When  shock  waves  that  are  not  jilane  (two-dimensional 
flow)  or  conical  (flow  with  axial  symmetry)  ocmir  in  uniform 
flow,  the  variation  of  entropy  across  the  shock  is  not  constant 
and  the  flow  hehiud  the  shock  is  no  longer  isentropic  and 
becomes  rotational.  If  the  variation  of  entropy  is  small, 
the  efleet  of  rotation  of  tin  flow  is  not  imjiortant  for  de¬ 
termining  the  pressure  distribution  along  a  body  and  the 
theory  of  potential  flow  gives  correct  I'csulls.  If  the  shock 
wave,  is  strong  and  has  large  curvature,  however,  th<‘  effect 


of  the  rotation  becomes  important  and  tlu'  flow  must  be 
considered  ns  rotational. 

The  method  of  characteristics  can  he  extended  to  apjily  to 
rotational  flow  if,  in  place  of  the  potential  function  for  the 
differential  equation  of  motion,  the  stream  function  con¬ 
sidered  by  Crocco  (reference  8)  is  used.  With  the  characti'r- 
istics  method  for  rotational  flow  a  more  exact  determination 
can  be  made  of  the,  shape  of  the  shock  wave  and  the  distribu¬ 
tion  of  velocity  and  pressure  for  phenomena  in  which  the 
effect  of  rotation  is  important,  as  in  the  internal  flow  through 
supersonic  inlets.  The  proeedure  of  numerical  calculation  is 
similar  and  not  much  more  complicated  than  that  used  foi’ 
the  case  of  potential  flow  with  axial  symmetry. 
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SYMBOLS 


pressure 

density 

cntro])y,  mechanical  units 

velocity 

Mach  number 

limiting  velocity  corres])onding  to  adiabatic  expansion 
to  zero  pressure 
speed  of  sound 


ratio  of  velocity  to  limiting  velocity 


x-component  of  relative  velocity 
7/-component  of  relative  velocity 
Cartesian  coordinates 

Macli  angle  ^arc  sin 

angle  between  velocity  f'and  x-axis,  radians 
angle  between  tangent  to  shock  and  direction  of 
velocity  of  flow  in  front  of  shock 
deviation  of  direction  of  velocity  across  shock  wave 
angle  of  polar  coordinate  in  conical  field 
potential  function 

stream  function  for  rotational  flow  (see  equation  (11)) 
ratio  of  specific  himt  at  constant  pressure  and  constant 
voltimc 
gas  constant, 
normal  to  streamline 
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.  _  sin  j3  tan  0  sin  ip 
”  cos  {<p+^) 

sin  8  tan  <3  sin  (f 
^  ~  cos  {<p—^) 

H,  L,  K,  and  N  defined  by  equation  (14) 

il  =  XA  —  XB 
io=Xc—XA 

,  Wa 


^  _  sin  ^A 
cos  {^a  +  <Pa) 

■_  sin 

^  cos  {(Pb  —  Pb) 
r=tan  /3^+A  tan  /3b 

_ L 

7^^  «+<;/ 

Subscripts: 

0  chamber  condition  (zero-velocity  adiabatic  transfor¬ 
mation) 

A  points  of  first  family 

B  points  of  second  family 

C  quantities  in  the  points  calculated  from  A  and  B 
X  derivative  with  resiicct  to  x 

y  derivative  with  resiicct  to  y 

a  ahead  of  shock 

6  behind  shock 

ipCn  value  corresponding  to  value  of  <p  at  point  C„ 

Cn  at  point  C„ 

CHARACTERISTICS  METHOD  FOR  SUPERSONIC  POTENTIAL 
FLOW  WITH  AXIAL  SYMMETRY 

Tlie  difTercntial  eciuation  for  potential  flow  with  a.xial 
symmetry  (reference  4)  is 


V  aVc)J--  V  a\)hy- 


a-  bxdy  ^y  ^  ^ 


In  supersonic  flow  some  lines  can  be  individuated  (charac¬ 
teristic  lines)  that  divide  the  flow  into  two  regions  for  wdiich 
the  values  of  <f>,  4>r.  and  <f>„  along  the  line  are  difrerenl.  For 
every  point  of  the  flow  two  characteristic  lines  can  be  <l('- 
tenninecr;  every  line  is  inclined  at  the  Mach  angle  with 
respect  to  the  direction  of  the  velocity  at  the  point,  aiul 
therefore  the  characteristic  lines  can  be  divided  into  two 
families  on  the  basis  of  the  sign  of  the  angle  of  the  characteris¬ 
tic  line  with  resju'ct  to  the  direction  of  the  velocity,  A 
family  that  is  usually  called  the  first  family  is  defined  by 
tiie  <‘()uati(  M 

V''=tan  (d+v^)  (2) 


and  the  otlu'r  family  (second  family)  is  defined  by 
(hi 

V—/3) 


The  variation  of  the  quantities  that  define  the  velocity 
(tp  and  F)  along  a  characteristic  line  is  given  by  the  following 
equations  from  reference  5: 

For  the  first  family, 

dp  tan  ^—1^=0  (4) 

and  for  the  second  family, 

+(Ip  tan  /3  — 77!  (5) 

wfierc  /  and  m  are  trigonometric  expressions  defined  as 

I  _  sin  8  sin  p  tan  g 
“  cos  (d+v) 

_sin  8  sin  p  tan  8 
cos  (ip-d) 

If  the  direction  of  the  velocity  and  the  Mach  number  at 
two  points  near  each  other  (points  A  and  B  in  fig.  1)  are 
known,  the  direction  of  the  velocity  and  the  Mach  number 
at  a  point  C  given  by  the  intersection  of  the  two  character¬ 
istic  lines  of  different  family  starting  from  points  A  and  B 
can  be  determined.  Because  the  distances  BC  and  J^C  are 
small,  all  the  coefficients  of  equations  (4)  and  (o)  can  be 
considered  constant  and  coincident  to  the  corresponding 
values  at  points  B  and  A.  and  the  tangents  to  the  character¬ 
istic  liiK's  at  points  A  and  B  can  bo  substituted  for  the  char¬ 
acteristic  lines  from  point  A  to  point  C  and  from  point  B 
to  point  0.  In  this  case  the  liiu's  AC  and  BC  are  straight 
lines.  The  line  BC  is  inclined  at  an  angle  p  —  8  with  respect 
to  tliex-axis  and  the  line  AC  is  inclined  at  an  angh^  ‘P-\-8  with 


X-  X 


KniCKK  I.— (Jpoinuirifal  corisirticiion  for  ik'trrniiiiinp  point  C  by  ilu'  chiirucli'riitiics  molhotl. 
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respect  to  the  x-axis.  Equation  (4)  can  be  applied  for  the 
line  AC  where 

dx=Xc—x,\ 

y=yA 

v=v, 

/=/.., 


and  equation  (5)  can  he  applied  for  the  line  BC  where 


If  a  stream  function  4'  is  assumed  to  he  defined  by  the 
following  equations  from  reference  8: 


the  equation 


4^,=  -yr{\-W^)y-' 


(11) 


(12) 


dx^Xc—fn 


y=yii 

r=r„ 

m  —  7n„ 


For  ])rae,tical  use.  etpiations  (4)  and  (5)  are  combined  and 
transformed  by  means  of  the  ratio  IF  of  the  velocity  1'  to 
the  limiting  velocity  V,.  This  ratio  is  defined  as 


=  (7) 


and  the  following  equations  are  obtained; 


is  a  function  of  only  the  stream  function  4'  (reference  8); 
and  therefore  from  the  equations  of  state,  continuity,  energy, 
and  steady  motion  the  following  equation  can  be  obtained; 

(j^-i)f{4d=o  (13) 

Equation  (13)  is  a  Monge-Ampere  equation,  and  if 


11=  \ 


u- 

c- 


K=- 


uv 

? 


dipA  (^tan  luu  0a^=1-||-^— (ys  — s^a)  tan 

(^c-a-B)-(a-c-x.)  ^  ^  (8) 

^^^  =  tan^W.A  +  ^-2^^  (9) 

Wc=WA  +  d\VA 

<f‘c=‘PA-‘rdip,, 

With  the  method  of  charaeteristies  (reference  4)  it  can 
be  shown  that,  if  a  deviation  of  a  streamline  which  wets  the 
body  occurs,  the  phenomena  on  the  corner  are  regulate<l  by 
the  same  laws  that  regulate  the  two-dimensional  flow; 
therefore,  the  tangents  of  every  ehai-acteristie  line  starting 
from  the  corner  are  Icnown.  If  the  initial  flow  conditions 
are  known,  the  step-by-stej)  calculation  of  all  the  ])hvsieal 
properties  of  the  flow  in  the  entire  field  is  jjermitted,  partic¬ 
ularly  the  calculation  of  the  sha])e  of  the  shock  wave,  and 
the  pressure  distribution  along  any  body  of  revolution 
with  axial  symmetrical  flow  in  cases  in  tvhich  the.  hypothesis 
of  potential  flow'  is  correct. 

characteristics  method  for  supersonic  rotational 
flow  with  axial  symmetry 

Supersonic  perfect  (low  is  rotational  when  the  flow  is 
l)reeeded  by  a  shock  wave  and  wlu'n  the  variation  of  entrojiy 
across  the  shock  wave  changes  from  point  to  point  behind 
the  shock.  In  this  ease  the  transformation  of  the.  fluid 
along  every  streamlim'  is  isentroi)ic  until  another  shock  wave 
occurs  in  the  fluid  (I'eferencc  8). 


(14) 

two  characteristic  families  with  tlic  following  equations  can 
be  obtained:  For  the  first  family, 


dy  K  L 

(15) 

dx~ii  yw  H 

(10) 

and  for  the  second  family. 

dy  K  IK~ L 

dx  iryn-  h 

(17) 

(18) 

If  M  is  the  normal  to  the  streamline,  equation  (11)  yields 

2 

4'^+4'^=f^V^  (l-lF-)^-> 


and 


curl  VxV  1  ,  1  ds 

- - ^^ylidn 


because  s  is  constant  along  every  streamline;  therefore 
7-1  1  ds  1 


.m= 


2y  R  dn 

'  t/ii  (1  — in)>-i 


(19) 

(20) 

(21) 
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and  the  following  e.xpressions  can  be  obtained: 

dy 


dx 


(  1  -  1I  ’=)  7-1  !• 


7  —  i 


—  7-f  2 

■lf=)7-i  W, 


dy, 


h.=yuA^-'^y^-)  7-i+^^«(i-W'^)7-. 


2yuW 
y—  1 


(22) 


When  equations  (21)  and  (22)  are  substituted  in  equations 
(lo)  and  (17)  and  (10)  and  (18)  and  the  Mach  angle  and  the 
veloeityiai'e  expressed  in  polar  coordinates,  etjuations  (If)) 
and  (17)  become 


tan  (/3+(p) 


(first  family) 


dy 

dx 


=  tun  {ip—fi)  (second  family) 


(22) 

(24) 


and  equations  (10)  and  (18)  become 


*  o  1  dx  ,  dx  dfi  sin“  n  t  e  -i  x  xo-x 

-  tan  ^  ^=0  (first  family)  (2o) 


</Tr  ,  .  .  ,  dx  dx  f/.' 

-jj. +tan/3(/9— - ;  7n  — 


sin^ 


2/  diiyRcos(ip—(i) 


=  0  (second  family)  (2C) 


Equations  (22)  and  (24)  are  identical  to  ecpiations  for 
potential  flow  (2)  and  (3),  and  equations  (25)  and  (20)  are 
similar  to  equations  for  potential  flow  (4)  and  (5),  differing 
only  bj'  the  terms  that  contain  dn.  Equations  (23)  to  (2C) 
permit  a  step-by-stej)  calculation  of  the  cntrojiy,  intensity 
of  velocity,  and  direction  of  the  velocity  if  the  initial  and 
boundary  conditions  are  knoiVn.  If  all  physical  properties 
are  known  for  two  points  A  and  B  and  if  the  two  points  are 
close  to  each  other,  the  tangents  to  the  characteristic  lines 
at  the  points  A  and  B  can  be  substituted  for  thi'  characteristic 
lines  with  close,  ajipi'oximation.  In  this  way  a  point  C  can  be 
determined  as  the  intersection  of  the  second  characteristic 
line  of  point  B  with  the  first  characteristic  line  of  point  A 
(fig.  1),  because  p  and  t3  are  Icnown  for  the  ])oints  A  and  B. 

For  the  characteristic  lines  of  the  first  family,  e<|uation 
(25)  gives  the  variation  of  ip  and  IF  from  jioint  A  to  iioint  Cf, 
and  all  the  coefficients  are  known  and  corresiiond  to  the 

coefficients  for  point  A :  only  the  term  is  unknown.  From 

equation  (21))  the  variation  of  ip  and  IF  from  point  B  to  jioint 

C  can  be  detei'mined,  and  all  other  terms.  exce|)t  ,  ai'c 

known  and  eipial  to  the  values  for  iioint  B.  The  term 

can  be  determined  from  the  value  of  the  entropy  for  points 
A  and  B. 

From  figure  1,  in  eiMiation  (25)  the  term  can  be  writ- 

di( 

ten  as 


dy _ t .Xf'  )  cos  (d.iTsc.i) 

dn  (.Cf.  — r., )  sin  d.-, 


(27) 


and  in  equation  (26)  can  be  written  as 

f/x  ^  ( Ac ^  '‘'c )^  cos  ( ipa  —  ^n) 
dn  (xc—xii'i  sin  d« 


(28) 


If  the  points  A  and  B  are  close  to  each  other  and  the  varia¬ 
tion  of  entropy  is  gradual,  equation  (27)  can  be  written  as 


ds 

dn 


d.y 


sin  d.i  I  ,  X  sin 
cos  (d.„  +  "  cos  (,,„-/3«) 


(29) 


and  can  be  considered  equal  in  equations  (25)  and  (26). 

In  this  case  is  known;  therefore  p  and  11  can  be  calculated 

for  point  C.  For  jiractical  calculations  equations  (25)  and 
(26)  should  be  transformed  into  two  eipiations  (equations  (31) 
and  (32))  each  of  which  contains  only  one  of  tlu'  unknown 
terms  dw  and  dp. 

For  sim])licity,  let 

h  =  T..|  —  Xb 


t(,=  Jc-F4 

h 


.'7=1  +  . 

«0 

i^‘P=‘P.l  —  <Pl! 

Bh 


h  = 


Wb 


sin  (3.1 


/■= 


cos  (^.I  +  I^.l) 
sin  ^B 


(30) 


cos  (pu  —  ^n) 
r=tan  tan  d.-i 

.1 

yJi  e+gf 

Then  the  following  equations  can  bo  obtained: 

rdp^=  1  — /(  — Aptan P„+A{yf  siir  d«  +  /ic sin-/3.i) 

I  inW  iig  _ 

?/»  y.i 
(/IF,, 


=tan  /3,,(/sp,,  — .lrsiir/3,,  +  7d' 

1 1 .1 


Pr=  Pa  -ddpA 

1Fc=1F,,  +  (/)F,. 

sur,3c=^T'(n',-l) 

In  equation  (31)  the  terms  and 

,'//<  ,'/.i 

portant  as  y  approaches  z('ro,  near  tin'  axis  of  the  body; 
therefore,  in  this  region  the  distance  between  the  points  con- 
sid(‘red  must  lie  reduci'd  to  obtain  sufficient  accuraejx 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

become  very  iin- 
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CHARACTERISTICS  METHOD  FOR  TWO-DIMENSIONAL 
SUPERSONIC  ROTATIONAL  FLOW 


In  two-dimensional  flow  an  eqmition  similar  to  equation 
(13)  can  be  obtained  if  a  special  stream  function  defined  by 
the  following  equation  is  assumed  (reference  8): 


In  this  case  the  equation  of  motion  (equation  (13))becomes 
(reference  8) 

( ^  “  cO  ~  i/'.. + ( 1  -  ■  2)  lA.. 

=  0  (38) 


Equation  (38),  like  equation  (13),  is  a  Monge-Ampere 
equation  and  permits  the  determination  of  two  equations  that 
define  the  characteristic  lines  and  two  equations  that  give 
the  variation  of  the  velocity  along  the.  characteristic  lines. 
The  equation  with  transformations  analogous  to  the  case  of 
three-dimensional  flow  can  be  written  in  the  following  form; 


^^=tun  (0+if)  (first  family) 

{ill 

^^=-tan  {(fi—0)  (second  family) 
dW  ^  dx  (Is  sin^  ^  ^  /c  ^  t  -i  i 

dW,  ^  „  j  dx  (Is  sin^  0  .  ,  j  r  -i  n 

IF  ^  ^‘^-yR  dn 


(39) 

(40) 

(41) 

(42) 


Equations  for  ip  and  dW  similar  to  ecpialions  (31)  and  (32) 
can  be  obtained  from  equations  (39)  to  (42)  by  using  equation 
(I  f{ 

(29)  for  the  expression  as  follows; 


r  dip,,=  1  —  h  —  Aip  tan  0„+A{(jf  sin-  0n+he  sin"  0,,)  (43) 


=  tan  0A  (lip—Ae  sin-  0,,  (44) 

Sr  =  Sj  +  AeyR  (45) 

‘Pr=  ‘PA~i~d>PA  (40) 

ire.=  ir,+(/ny  (47) 

sin-  0r-=^  .A  '  (]jy2-  l)  (4-3) 


DETERMINATION  OF  SHAPE  OF  SHOCK  AND  PRESSURE 
DISTRIBUTION  ALONG  A  BODY 

The  physical  properties  of  supei-sonic  flow  past  a  body  of 
revolution  in  axially  symmetrical  flow  can  be  determined 
step  by  step  by  the  use  of  equations  (31)  to  (3()).  The  sys- 


I-'KJfUE  2.— l‘mctic:il  sysU'TUof  culculalinp  Ihu  flow  fluhl  for  n  pointod-nosu  fxxly  of  revolution. 


tern  of  calculation  is  similar  to  that  for  potential  flow.  If  the 
body  begins  with  a  point  (reference  5).  a  cone  tangent  to  the 
body  (fig.  2)  can  be  substituted  for  the  front  part  of  tlie  body. 
If  point  A  is  the  point  at  which  the  body  can  no  longer  be 
considered  coincident  witli  the  cone,  the  velocity  at  imint  A 
is  known  from  the  cone  cahndations  (referenei's  9  and  10); 
therefore,  the  shaju'  of  tin'  characteristic  line  AH  of  the  first 
family  can  be  designed,  beeaiise  </>  and  0  corresponding  to 
different  angles  6  are  known  from  tin-  cone  calctilations. 
-\t  point  A  the  body  turns  through  an  angle  Ap  and  the  flow 
undergoes  a  transformation  that  is  determinable  by  the  laws 
of  two-dimensional  flow.  The  velocity  and  direction  of  the 
flow  after  the  deviation  Aip,  therefore,  and  the  tangiuits  to 
the  new  characteristic  lines  of  the  first  family  at  point  A  can 
be  designeil.  At  a  point  B,  near  point  A,  the  intensity  and 
direction  of  the  velocity  are  known;  consequently,  tlu^  tan¬ 
gent  to  the  characteristic  line  of  the  second  family  can  be 
designed  at  point  B,and  the  point  C,  can  lx-  located.  At 
point  C,  the  physical  properties  can  be  calculated  with  the 
equations  of  potential  flow  (equations  (8)  and  (SO )  because 
tile  shock  in  front  of  the  body  is  coniiail.  With  tin-  same 
system  the  point  C„  is  determined  from  the  jioint  B„  on 
the  shock  wave.  In  order  to  determine  the  flow  of  the  jioint 
C„4,  on  the  shock  wave,  the  equation  across  the  shock  and 
the  equation  for  the  cliaracteristie  lines  of  the  first  family 
must  be  used.  If  t  is  the  angle  Ix-tween  the  tangent  to  the 
shock  and  the  direction  of  velocity  of  flow  in  front  of  the 
shock,  6  the  deviation  of  the  diri-ction  of  velocity  across  the 
shock  wave,  and  the  subscrijits  a  and  b  denote  the  conditions 
ahead  of  and  behind  the  shock,  rcsjx-ctively,  the  shock  (-([na¬ 
tions  can  be  written  in  the  following  form  (reference  11); 


tail  t 


tan  (e— 6)  7H-  1  L^\/h‘  sin-  (t  — 5)~^  2  J 

tan  6  \  2  .1/,,-  sin-  e— 1  /  ' 


(49) 

(50) 
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SO 

(377^111^+ V)  (^2) 

(ios.  f^+7  loll,  "”)  (5:5) 

7  ^  Pa  Pft/ 

1  ,  1  2  1 

11-='+^-!  Mu  ^^4) 

Equations  (51)  to  (.53)  show  that  tlio  ficnpral  ponditioo  for 
obtaininft  potential  flow  beliiiul  tlu'  shock  when  the  flow  in 
front  of  the  shock  is  potential  is  that  M„  sin  «  must  be  con¬ 
stant  .  Since  tbe  value  of  M„  in  front  of  the  shock  is  known 

from  equations  (49)  to  (54),  the  values  of  -''’*.7''"^  and 

«o  <16 

can  be  determined  as  functions  of  6  and  therefore,  the  equa¬ 
tion  for  the  characteristic  lines  of  the  first  family  (equation 
(25))  becomes 


(/ir  ,1^ 

ll'c,  '^~(I6  U’c 


I  u  n  uip  ,  ax 

—  If <h-„-  U- 


+  — yp  +  rfs=]=0  (55) 

where  11’,^^  is  the  velocity  behind  the  shock,  s,_,  is  the 

ontrop}',  and  are  coofliciiuits  that  corres|)ond  to 

a  deviation  across  the  shock  for  which  the  velocity  assumes 
the  direction  of  (pc,,,  H  c„  is  the  velocity,  and  is  the  entropy 
at  the  point  C„.  When  the  value  of  </(p  has  been  detcrinined, 
the  deviation  across  the  shock  6  and  the  correspondiiif;  values 
of  W  and  «  can  be  determined  at  point  C„+|, 

In  order  to  determine  the  velocity  on  the  body  at  a  point 
A,,  equation  (2G)  is  used.  At , point  A,  the  value  of  p  is 
known  because  the  direction  of  the  flow  is  tangent  to  the 

body  and,  therefore,  dip  is  known;  is  also  known  because 

in  equation  (27)  the  value  of  Sc  corresponds  to  the  value  of 
s  at  the  point  Aj  and  is  equal  to  the  value  at  point  A,  which 
is  known  from  the  eone  calculation.  From  the  value  of  IF, 
th('  pressure  on  the  body  relative  to  the  pressun'  p,,'  (pressure 
for  zero  velocity  from  isetitropic  transformation  from  the 
conditions  behind  the  shock)  can  b('  deti'rmined  as 


,  =  (1  — If'-)  7-1 


On  the  lip  of  the  nose  of  an  ojx'ii-no.se  body  a  shock  wave 
that  is  a  two-dimensional  shock  occurs;  therefore,  the  tan¬ 
gent  to  the  shock  on  the  lij)  is  known  and  the  jiressnre  and 
V(>locity  behind  the  shock  ai'e  also  known  (eciuations  (49) 
to  (.54)).  In  figure  ;5  the  line  AB  and  the  velocity  l)chind 
tl;  '  shock  at  point  15  are  known.  With  equation  (20)  the 
^ t4o(‘ity  at  point  C  can  be  calcnlateil  with  tlu'  svsti’in  just 
described,  and  from  th(>  point  B  the  point  D  along  the  charac- 
teii.'.tic  line  of  the  first  fainilv'  can  be  calculated  by  equation 


(55).  From  jioint  D  the  point  E  can  be  calculated  and  a 
point  F  can  be  interpolati'd,  which  permits  the  determination 
of  points  G  and  II.  Because  the  curvature  of  the  shock 
near  the  lip  is  large,  a  point  C  that  is  very  near  point  A  must 
lie  selected  so  that  correct  interpolated  values  mav’  be  ob¬ 
tained.  The  point  E  can  be  recalculated  from  point  G,  and 
the  points  I  and  M  can  be  recalculated  from  point  F  in  order 
to  have  a  second  apjiro.vimation. 

For  the  case  of  two-dimensional  flow  the  procedure  is  the 
same  as  that  for  three-dimensional  flow;  equations  (43)  to 
(48)  are  used  instead  of  the  corresponding  e(|uations  (e(|ua- 
tions  (31)  to  (30)). 

E.XAMPLES  OF  APPLICATION  OF  CHARACTERISTICS 
METHOD  FOR  ROTATIONAL  FLOW 

nETERMlNATION  OK  SHOCK  SHAPE  AND  PRF.SSl'RE  DISTRlIitlTION 
ALONG  THE  EXTERNAL  SURFACE  OF  A  SLENDER  OPEN-NOHE  BODY  OF 
REVOLUTION 

Theoretical  and  experimental  calculations  were  made  to 
di'termine  the  shock  shaiie  and  pressure  distribution  along 
the  external  surface  of  a  slender  open-nose  body  of  revolu¬ 
tion.  The  body  considered  is  the  sami'  bodv  for  which  cal¬ 
culations  of  the  external  pressure  distribution  were  made  by 
Brown  and  Barker  of  the  Langlev  Memorial  Aeronautical 
Lalioratory  by  use  of  tbe  small-disturbance  theory.  Tlu' 
calculations  were  made  for  a  free-stream  Mach  number  of 
1..525,  for  which  a  schlieren  photogra])!)  taken  during  testa 
was  available  for  eonqiarison  of  tlu>  calculated  and  test 
results. 

In  order  to  determini'  the  importance  of  rotation  of  the 
flow,  calculations  for  the  front  part  of  the  body  were  also 
made  with  the  potential-flow  characteristic  equations.  The 
results  of  these  calculations  differed  only  slightly  from  those 
obtained  from  the  characteristics  method  for  rotational  flow 
because  the  curvature  of  the  shock  was  very  small.  The 
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FKjI'RE  4.— Practical  system  of  calculiUinp  the  external  flow  around  a  slender  o|M‘n*no.se  lH)dy 
of  revolution  for  1.525. 


calculations  wore  begun  with  the  determination  of  the  two- 
dimensional  shock  on  the  lip  of  the  nose.  The  practical 
system  of  the  calculations  is  shown  in  figure  4,  and  in  figure 
5  the  calculated  streamlines  and  shock-wave  shape  are  com¬ 
pared  with  the  shock-wave  shape  obtained  from  test  residts. 
In  figure  6  the  pressure  distribution  calculated  by  the  char¬ 
acteristics  method  is  compared  with  the  pressure  distribution 
determined  by  the  small-disturbance  theory.  The  small- 
disturbance  theory  undervalues  the  increase  in  pressure  that 
occurs  through  the  shock,  but  the  differences  in  the  results 
obtained  by  this  method  and  those  obtained  by  the  charac¬ 
teristics  method  are  small. 


o  Experimental  o 


Kir.UHE  5.— Calculiiiod  .slntamliiifs  and  shock  wave  for  a  .slender  oix'n-nosr  body  of  revolution 
for  A/ ■=  1.525,  showinp  exixjrimcntal  .shock-wave  shajHb 


I>ETERM1NAT10N  OF  SHOCK  SHAPE,  STREAMLINES.  AND  PRESSURE  DIS- 
TRIBUTION  ALONG  THE  INTERNAL  SURFACE  OF  A  SLENDER  OPEN-NOSE 
BODY  OF  REVOLUTION 

Three  slender  open-nose  bodies  with  different  conical  inlet 
.angles  are  considered  (figs.  7  to  12)  and  the  supersonic  part 
of  the  internal  flow  is  studied  for  a  free-stream  Mach  number 
of  1.0.  For  this  type  of  body  the  internal  shock  produced 
at  the  lip  of  the  inlet  has  a  very  large  curvature  and  the  effect 
of  rotation  is  therefore  very  imj)ortant.  The  calculations  are 
extended  to  the  region  in  which  the  Mach  number  is  1.0. 


Figi  reO,-  Pressure  (iistrihution  alonp  the  cxtornnl  surface  of  a  slender  oinm-mist' body  of  rovolulion  for  A/=1. 525.  (I)ntn  for  small-disturhaneo  theory  from  work  by  Urown  and  Parker 

of  the  Lanpley  Laboratory.) 
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Figvre  7.— Practical  system  of  calcuhuinp  the  internal  sui)ersonic-flow  (juantities  for  a 

slender  o{K‘n-nose  body  of  revolution  for  Kiouke  8.— Calculated  streamlines  and  shock  wave  for  internal  flow  in  a  slender  oiK'n-nose 

body  of  revolution  for 
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Figure  ll.— 'Practical  sysU'in  of  culculuiinR  iho  Internal  su}x‘rsonlc*flow  proix^rtlcs  for  a 
slender  o|K*n-no.se  body  of  revolution  for 

The  hypothesis  is  made  that  subsonic  boundary  conditions 
and  stability  considerations  permit  a  subsonic  flow  at  the 
end  of  the  supersonic  flow  such  as  results  from  the  calcula¬ 
tions.  The  results  of  the  calculations  show  that  at  the 
axis  of  the  inlet  a  normal  shock  always  occurs  and  that  the 
region  in  which  a  strong  shock  occurs  (with  subsonic  velocity 
behind  the  shock)  increases  in  extension  with  the  increase 
in  angle  of  the  internal  cone.  When  the  cone  angle  aj)- 
proaches  zero  and  the  shock  is  a  Alach  wave,  a  complete 
reflection  occurs  at  the  axis  of  the  inlet  and  the  extension  of 
the  strong  shock  is  zero. 

For  large  internal  cone  angles  (figs.  9  to  12)  the  shock  is 
a  simple  shock  that  liecomes  normal  in  the  central  part  of 
the  body  of  revolution.  After  the  shock  the  compression 
continues  but  the  characteristic  lines  cannot  form  an  enve¬ 
lope.  For  small  internal  cone  angles  compression  occurs 
gradually  and  an  envelope  of  Alach  lines  occurs.  The  shock 
therefore  reflects  from  the  central  part  and  another  shock  is 
generated.  The  form  of  the  reflected  shock  is  shown  in  figure 
8.  The  ratio  of  the  diameter  of  the  region  in  which  a  strong 
shock  occurs  to  the  diameter  of  the  inlet  as  a  function  of  the 
internal  cone  angle  is  shown  for  i\/=1.0  in  figure  i:5.  Tin* 
results  are  interesting  for  the  practical  design  of  supersonic 
inlets  of  slender  shape  hecause  they  show  that  for  large  cone 
angles  the  central  part  of  the  body  of  revolution,  in  which 
the  compression  is  not  v('ry  efficient,  is  large. 


Figure  J2.“-ralcul«tt‘d  str<*amliiifs  and  shock  wave  for  intornul  flow  in  a  slcndor  o|x*n*nose 
body  of  rovolullon  for 

DETERMINATION  OF  PHYSICAL  PROPERTIES  OF  THE  INTERNAL  FLOW 
THROUGH  AN  INLET  WITH  A  CENTRAL  BODY 

For  a  Alach  number  of  l.G  an  analysis  of  the  shape  of  an 
hdet  with  a  central  body  was  made  to  aid  in  obtaining  high 
efficiency.  Theoretically  a  supersonic  diffuser  with  or  with¬ 
out  a  central  body  and  having  no  shock  losses  or  shock  drag 
can  be  obtained  (reference  12);  but  for  practical  use  it  is 
convenient  to  accejjt  small  sliock  drag  in  order  to  avoid  large 
friction  drug.  The  inlet  considered  (fig.  14)  has  a  10°  central 
cone.  The  deviation  aeross  the  conical  shock  is  48'.  The 
shock  is  reflected  by  a  cylinder  that  forms  the  extc'rnal  part 
of  the  inlet.  The  reflected  shock  produces  rotational  flow 
behind  the  shock  and  the.  deviation  across  this  shock  on  the 
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Figube  13.-“Kaiio  of  the  diameter  of  the  reirion  in  which  strong  sliock  occurs  tt>  the  diameter  of  the  inlet  ns  a  function  of  the  angle  of  internal  cone.  A/=1.0. 


I  louRf;  1-t.  rtiictieal  system  of  caleulnt ing  the  shape  of  the  central  body,  of  the  streanilini's.  and  of  shock-wave  shape  for  a  supersonic  inlet  at  -\/=l.G. 
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cyliiuiiT  is  1®  and  on  the  contral  liody,  2°  43'.  If  the  con- 
trul  body  bcliind  tlic  shock  lias  tlic  same  direction  as  the 
velocity,  the  shock  will  not  be  reflected  and  isentropie 
compression  can  be  obtained  behind  the  shock  (fi^.  14). 
The  design  of  the  central  body,  therefore,  was  deter¬ 
mined  from  the  calculation  of  the  corresponding  streamline. 
The  variation  of  the  velocity  along  the  external  cylinder  and 
the  value  of  the  exit  velocity  were  fixed,  and  from  this  con¬ 
dition  and  the  condition  dependent  on  the  shock,  the  ve¬ 
locity  at  every  point  was  calculated.  In  order  to  avoid 
errors  stream  tubes  were  designed  that  permitted,  on  the 
basis  of  the  ratio  of  the  area  in  the  region  of  uniform  velocity 
to  the  area  at  the  end  of  the  stream  tube,  a  check  on  the 
precision  of  the  numerical  calculations.  Th(‘  -Mach  num¬ 
ber  in  the  minimum  section  of  the  inlet  was  fixed  at  a  value 
larger  than  1.0  so  that  distui'bances  from  the  subsonic  part 
of  the  How  would  not  caus('  instability.  The  value  chosen 
was  1 .00. 

If  an  effuser  is  connected  with  the  diffuser,  a  body  of 
revolution  with  very  low  shock  drug  can  be  obtained  (fig. 

15  i.  The  only  pressun'  losses  are  the  losses  across  the  two  ' 
shocks,  which  are  v<'ry  small:  but  for  piactical  a|)i)licat.ions  | 
tb('  friction  losses  are  largei-  tban  for  tlu'  internal  bodv  alone.  I 
A  balance  of  tbe  i)i'essure  losses  and  friction  losses  must 
tluT<'foi'e  be  made  in  ordei-  to  examine  tbe  |)ossibility  of  i 
jiraetical  use  of  this  arrangement, 

CONCLUSIONS  I 

A  system  foi-  caleulaling  tbe  jibysical  ]n'opei-ties  of  super-  i 
sonic  rotational  flow  with  axial  symmetry  and  suiiersonic  ' 
rotational  flow  in  a  two-dimensional  field  was  determined  . 


by  use  of  the  characteristics  method.  Practical  use  of  the 
system  is  based  on  a  stc])-by-step  procedure,  which  requires 
long  numerical  calculations;  but  the  calculations  for  three- 
dimensional  flow  arc  of  the  same  type  ns  for  jtotential  flow 
and,  th(>rcfore,  can  be  used  for  the  practital  jiroblems  in 
which  rotation  is  important.  Some  applications  were  made 
to  determine  the  external  and  internal  flow  on  bodies  of 
revolution  with  axial  symmetr}',  and  the  following  con¬ 
clusions  were  indicated: 

].  The  effect  of  rotation  is  not  very  important  if  the 
variation  of  entropy  is  small  but  is  important  in  the  study 
of  internal  flow,  for  which  the  variation  of  entropy  is  usually^ 
large. 

2.  When  the  inlet  is  a  circular  conical  channel,  a  shock  is 
produced  at  the  lip  of  the  inlet  that  becomes  stronger  when 
the  shock  approaches  the  axis  of  the  mlet  and  becomes  a 
normal  shock  at  the  axis.  The  region  in  which  a  strong 
shock  occurs  (with  subsonic  velocity  behind  the  shock)  in¬ 
creases  with  the  inci'easi'  of  the  angle  of  intenial  cone. 

3.  If  an  inlet  with  a  central  body  is  considered,  the  method 
of  characteristics  permits  the  design  of  an  internal-channel 
sha])e  that,  theoretically,  lesults  in  very  eflieieiil  recom- 
pi'ession  in  the  inlet :  and,  if  an  effuser  is  conneeled  with  the 
diffus(“r,  a  body  of  revolution  with  very  small  sboek-wave 
drag  can  be  detei'inined. 


Langley  Memouial  Aekonaftical  Lahouatoky, 
National  Advisory  Committee  for  Aero.naftics, 
Langley  Fiei.d,  \'a.,  April  29,  1949. 


12 


RKPOHT  NO.  S  ll  —  XATIOXAL  ADVISOIiA’  COMMITTER  I'OIi  AEHONACTICS 


KEFERENCES  ! 

1.  Praiidtl,  L.,  and  liusciiiann.  A.:  Xalinrun^svcTfahron  ziir  znicli- 
nerischon  Krinlttlung  von  obcnon  Sirdinnii'joii  mit  Dlxirschall-  I 
^eschwindi^keil.  Stodola  Festschrift,  Zliricii  and  Loipzij;,  1929.  | 
*2.  Ihiseniann,  A.:  Ciasdynanjik.  Ilandb.  d.  Fx))('rini(‘nl.alpi)ys.,  Bd.  i 
IV,  1.  Tcil,  Akad.  Wriatrscesellschaft  in.  1).  IF  (Leipzig),  1931,  | 
pj).  343-4G0.  i 

3.  Frankl,  I‘\  F:  Sujiersonic  Flows  with  Axial  Svinmotry.  Izvcstiu 

Akadeinii  RKK.V,  vol.  1,  1934. 

4.  Ferrari,  C.:  Cainpo  aorodinainico  a  velociia  iperaoustica  attorno 

a  un  solido  di  rivoluzione  a  prora  acuminata.  F'.-Verotecnica, 
vol.  XVI,  fasc.  2,  Fot).  HISO,  pj),  121-i:«).  j 

5.  Ferrari,  The  Determination  <tf  the  Projectile  of  Minimum 

^\'ave-Ucsistallce.  IF  T.  P.  'I’ranslat ion  No.  1189,  British 
Ministry  of  .\ircraf(  Production.  (Pan  II  from  Atti  IF  Acead. 
Sci.  3'orino.  vol.  7o,  Xov.-Dec.  1939,  {>]).  (5 1-99.)  | 

6.  (iuderley,  (F:  Die  (lharakteristikeiimethode  fur  olxme  und  j 

achsensymmelrisch(“  Dberschallstromunpen.  Jahrb.  1940  der  | 


deutscheti  Lufifahrlforschunji,  IF  Oldenboiira  (MunicliF 
pp.  I  522  -  I  535. 

7.  Sauer,  P.:  Charakterislikenvcrfahren  fur  raiimlielie  achsen.syin- 

meirische  tlberschallstromungoii.  Fursrhungsix'richl  Xr.  1299, 
I)eutsch(‘  Fuftfahrlforschung  f.Vachen),  1940. 

8.  Crocco.  L. ;  Una  mtova  funzioiie  di  correme  jter  lo  studio  del  moto 

rotazionale  dei  aas.  Rend,  della  R.  Accad.  dei  Fincei.  ^'oi. 
XXIII,  -^er.  9,  fasc.  2,  Feb.  1939. 

9.  Busematin,  Driicke  auf  kegelformige  Sf)iizen  Ixu  Beweguiig 

mit  'Cberschailae.-chwindigkeit.  Z.  f.  a.  M.  ^F,  Bd.  9,  Heft  9, 
Dec.  1929,  pp.  499-498. 

10.  Taylor,  (1.  F.  aud  Maccoll,  .1.  \V.:  The  .\ir  Pn'ssun!  on  a  Cone 
Moving  at  High  Sf)eeds.  Proc*.  Roy.  Soc.  (Londoti),  s(‘r. 
vol.  139,  no.  838.  Feb.  1,  1933.  j)p.  27S-31F 
IF  (j‘occo,  Luigi:  Sinaoiarit:).  della  corrente  gassosa  iperaeust lea 
neiriniurno  di  una  jtrora  a  dicdro.  F’.AerotecJiica,  vol.  X\’II, 
fa.sc.  9.  June  1937,  pjn  519-534. 

12.  Ferrari,  ('amj)i  di  cornmie  ipersonora  attorno  a  soiidi  di 
rivoluzione.  L'Aerotecnica,  vol.  XN'II,  fasc.  9,  Juik*  1937, 
])p.  507-518. 


U.  S.  GOVERNMENT  PRINTING  OFFICE:  1949 


Y 


X 


Positive  directions  of  axes  and  angles  (forces  and  moments)  are  shown  by  arrows 


Axis 

Force 
(parallel 
to  axis) 
symbol 

Moment  about  axis 

Angle 

Velocities 

Designation 

Sym¬ 

bol 

Designation 

Sym¬ 

bol 

Positive 

direction 

Designa¬ 

tion 

Sym¬ 

bol 

Linear 
(compo¬ 
nent  along 
axis) 

Angular 

X 

X 

L 

Y - tZ 

Roll  ... 

u 

P 

Lateral..^. 

Y 

Y 

M 

Z - tX 

Pitch. _ 

9 

V 

Q 

z 

z 

N 

X — ►y 

Yaw _ 

4' 

tD 

r 

Absolute  coefficients  of  moment 

o-J^  c-M. 

^~^cS  '^”~qbS 

(rollmg)  (pitchmg)  (yawing) 


Angle  of  set  of  control  surface  (relative  to  neutral 
position),  5.  (Indicate  sirnface  by  proper  subscript.) 


4.  PROPELLER  SYMBOLS 


D  Diameter 

p  Geometric  pitch 

pjD  Pitch  ratio 

V'  Inflow  velocity 

V,  Slipstream  velocity 

T  Thrust,  absolute  coefficient  Ct— — Sn 

pvrD* 

Q  Torque,  absolute  coefficient 


P 

C. 

V 

n 

# 


Power,  absolute  coefficient  Cp— 


P 


Speed-power  coefficient = 


Efficiency 

Kevolutions  per  second,  rps 


Effective  helix  angle = tan 


6.  NUMERICAL  RELATIONS 

1  hp=76.04  kg-m/s=550  ft-lb/sec  1  lb=0.4536  kg 

1  metric  hor8epower=0.9863  hp  1  kg=2.2046  lb 

1  mph=0.4470  mps  1  mi=  1,609. 35  m=5,280  ft 

1  mp8=2.2369  mph  1  m=3.2808  ft 


